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D 

tin- Abstract. We consider irrotational dust spacetimes in the full non-linear regime which are 

■ "purely radiative" in the sense that the gravitational field satisfies the covariant transverse 

^SJ ' conditions div(_ff) — div(_E) — 0. Within this family we show that the Bianchi class A spatially 

homogeneous dust models are uniquely characterised by the condition that H is diagonal in the 

' shear-eigenframe. 

> 

o 

1. Introduction 

, Irrotational dust spacetimes are typically considered as models for the late universe or for 

I gravitational collapse. They are covariantly characterized by the dust four- velocity u"", energy 

O . density p, expansion 9, and shear aab, and by the free gravitational field, described by the electric 

and magnetic parts of the Weyl tensor Cabcd- 

U 

bX). _ c d cd e 

Eab = CacbdU^U — E^ab> Hab = -^^acdC^ beU^ — H^ab> (1) 



X 
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where Sabc = Vabcd'^J-'^ is the spatial projection of the spacetime permutation tensor r]abcd and 
^ i S<ab> = ha'^hj,'^S(^cd) ~ ^Scdh'^^hab is the projected, symmetric and trace-free part of Sab with 

hab = dab + Ua^b the Spatial projector into the comoving rest space, and gab the spacetime 
metric. Gravitational radiation is covariantly described by the nonlocal fields Eab, the tidal 
part of the curvature which generalizes the Newtonian tidal tensor, and by Hab, which has no 
Newtonian analogue [5j. As such, Hab may be considered as the true gravity wave tensor, since 
there is no gravitational radiation in Newtonian theory. However at least in the linear regime, 
as in electromagnetic theory, gravity waves are characterized by Hab and Eab, where both are 
divergence- free but neither is curl- free [2]. Here the covariant spatial divergence and curl for 
tensors are defined by [3]: 

{dwS)a = D'^Sab, cnvlSab = S^aiaD'S^/, 
where the spatial part Da of the covariant derivative is given by 

DaAi,...=ha'h'^---VcAd... 

For this reason spacetimes in which div(£') = = div{H) have been termed purely 
radiative p!]. Although, when considering the gravitational dynamics of matter in the full non- 
linear regime, these conditions give rise to two independent chains of rather severe integrability 
conditions, it was remarked in [T] that both are satisfied by the Bianchi class A solutions and 
hence purely radiative solutions at least can exist. We turn the argument around and show that. 



within the family of irrotational dust spacetimes for which div(£') = = div(H), the Bianchi 
class A spatially homogeneous dust models are uniquely characterised by the condition that H 
is diagonal in the common shear and E-eigenframe. The existence of the latter eigenframe is 
guaranteed by the Bianchi identity ([8]) below, 

where the spatial dual of the commutator of tensors, [S,T]a, is defined as [/S, Tj^ = EabcS^'^Tif^ . 



2. Basic variables and equations 

We follow the notations and conventions of the Ellis-MacCallum orthonormal tetrad 
formalism [4j. Herein, the normalized 4- velocity u*^ plays the role of the timelike basis vectorfield 
Cq of the tetrad {cq, ej^}. Greek indices take the values 1, 2 and 3 and refer to tetrad components 
with respect to e^. 

The basic variables in this formalism are 24 independent linear combinations of the Ricci 
rotation coefficients or, equivalently, of the commutator coefficients associated with the 
tetrad; these are the objects Uap and Oq defined by = epse'n'^"' + 5faj3 — S^^as, the components 
Via of the angular velocity of the triad e^ with respect to the local 'inertial compass' and the 
components iia, i^a and 9^^ of the spatial kinematic quantities (acceleration vector, vorticity 
vector and expansion tensor respectively). We also use the decomposition O^p = (Tap + 1/S65af3, 
with fjafe the trace- free shear tensor and 6 the expansion scalar. It turns out that there is a slight 
notational advantage in using instead of the variables riap and the connection coefficients n^, 
Qa and related to the former by 

na+l Q-1 = {ra + ga)/2, a„ = (r„ - qa)/2, Uaa = ria+i + na~i 

(these expressions have to be read modulo 3, so for example a = 3 gives ni2 = 2^3 + rs). 
For sake of simplicity we define and rria as the spatial gradients of 6 and p respectively: 

Za = daO, nia = daP- 

Together with the Jacobi identities and the (a/3) field equations [4], the basic equations in the 
formalism are the Ricci equations for Ua and the Bianchi equations [5], which for = Uq, = 
simplify to the following: 

Ricci equations 

dWaa - "^DaO = 0, (2) 

curluab - Hab = 0, (3) 
2 

(T<ab> + -^d(Tab + <Tc<aCrb>'' + Eab = 0, (4) 

e + ^e^ + aaba'^' + l{p + 3p) = 0, (5) 

Bianchi equations 

p = -{p+p)e, (6) 

divEa = [a,H]a + ^DaP, (7) 
diyHa = -[a,E]a, (8) 
E<ab> - cmlHab = -OEab + Sa^^aEby" - ^{p + p)crab, (9) 



H<ab> - curl£;„fc = -9Hab + 3ac<aHb>''. (10) 
As we are working in the {a, i?)-eigenframe we derive from (jl]) with a ^ (3 that, 

f^i(o-22 - 0-33) = ^^2 (0-33 - 0-11 ) = ^^3('7ii - 0-22) = 0. (11) 

As a consequence the angular velocity Via can be assumed to be 0: either all components 
automatically vanish, or, when e.g. the (62,63) plane is a shear eigen-plane, then 0.2 = ^^ = ^ 
and we can choose an extra rotation in this plane to make Jli = 0. The frame is then fixed up 
to rotations in the (62, e3)-plane for which the rotation angle 99 satisfies Oq^p = 0. 

From the diagonal components of ([3]) we obtain the following algebraic relations between 

-^aa^^cta and TIq- 

Hii = n3((T22 - 0-11 )+ f^2 (0-33 - 0-11 ), 

H22 = ni((T33 - 0-22) + n3((Tii - CT22), (12) 

-f^33 = mi^^ii - 0-33) + "-1(0-22 - 0-33) • 

From the definition of divi^^ and the assumption of div(£') = we have 

diEii = Ellin -2qi)- E22{ri +qi), 

52-E22 = -©22(^-2 - 2^2) - £^33(?'2 + Q'2), (13) 
53S33 = S33(r3 - 2^3) - ^ii(r3 + ^3), 



while the Bianchi equation ([7]) yields 



mi = 3-^23 (o"33 - 0-22), 

m2 = 3iJi3(crii - 0-33), (14) 

771-3 = 3i?i2(fT22 - O-ii). 

From ([2]) we 03-11 derive the spa-ticil derivatives of cTcxol 

diuii = + ri(o-ii - CT22) - gi(o-ii - 0-33), 

<920-22 = 1^2 + ^2 (-722 - 0-33) - 92(0-22 - 0-11 ), (15) 
•930-33 = 3^3 + ''3(0-33 - 0-11) - 93(0-33 - 0-22)- 

Note that we can assume that the shear is non-zero, as otherwise solutions would be conformally 
flat and we have a FRW universe. 

We now act with the commutator [c?o, da] on p to obtain the relation 

4 

doma = -ma{aaa + -0) - {p + p)Za., (16) 



after which the time evolution of (jl4p results in 
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3(aii - a22)doHi2 = {6a(^ - 6ai^ - QE22 - 3^33 + 2^22^ - 2aii6)Hi2 + {p + p)Z: 

3(^22 - a33)5oi/23 = (6^12 -6a|3-6£;33- 3^11 + 2(733^- 2(7220)^23 + (P + P)^1 (17) 

3(^33 - (Tii)doHvi = (6ai3 - Qali - QEu - 3^22 + ^auO - 2a^ze)Hi^ + [p + p)Z2 

Provided we assume that also H^/s is diagonal in the {a, £')-eigenframe, it follows directly 
from ([III) and (fT7|) that Za = and daP = 0, suggesting that the corresponding spacetimes 
might be spatially homogeneous. Furthermore, as [3o,9a] = —{(7aa + ^0)da we obtain from ([5]) 
and daO = that daO'abO'"'^ = and hence, using ([3]) and (fTSl) . 

n(o-ii - 0-22)^ - 91(0-11 - 0-33)^ = 

^2(0-22 - 0-33)^ - 92(0-22 - 0-11)^ = (18) 
^•3(0-11 - 0-33)2 - 93(0-33 - 0-22)^ = 0. 

We will treat the two cases of non-degenerate shear and degenerate shear now separately. 



3. Non degenerate shear 

Propagating the last three equations along the fluid flow we get 

gi(<Tii - Cr33)x = ^2(0-11 - 0-22)x = «3(0-22 - 0-33)X = (19) 

where x is defined as 
2 

X = -^[Ell{a22 - 0-33) + £'22(0'33 - C^ll) + -E'33(<7ll - Cr22)] - ((Tn - (T22)(c^22 " <733)(f'"33 " <7ll). 

If (?f + ^2 + 93 7^ it would follow from (fT9l) that X = 0) after which cJix = leads to the 
contradiction (7i(Tii((Tii — (T33) = 0. Therefore we conclude that = and, according to 
([TH]) . we find that = 0, such that = and n^^ is diagonal. It is easy to check that all 
spatial derivatives of the rotation coefficients now vanish and, as the frame is invariantly defined, 
we obtain the non-degenerate spatially homogenous models of Bianchi class A (i.e. types 
I,II,Vh,VIIo,VIIIJX). 

4. Degenerate shear 

From the degeneracy of the shear it is easily verified that Eap must be degenerate too; so we 
can assume 

Gap = diag(-2(T, fj, a) 
Eaf3 = diag{-2E,E,E). 

Assuming that Ha/3 is diagonal in the {a, ii^)-eigenframe, (jlSp implies 

ri - qi = rs = q2 = (20) 

while the (23) component of ([3]) shows that ri + gi = and hence ri = qi = 0. 

At this point the only variables having possibly non- vanishing spatial derivatives are 93, r2 
and ni, 71-2,^3. Eliminating E fromQ for a = (3 = 1 and a = /? = 2, we obtain the following 
equation 

= Hii{-2ni +n2- 2n^) + i?33(-4ni - 722 - 113), 
In combination with ()12p the latter leads to 

(ni +712 + n3)(n2 - 713) = 0. (21) 

We will discuss the three cases that follow from (|2T]l separately; i.e. n2 7^ n3, 712 = W3 7^ and 
^^2 = ^^3 = 0. Note that in the ensuing classification E is nowhere zero, due to the non-existence 
of so-called Anti-Newtonian universes ^ 

4.1. 712 - ri3 7^ 

In this case ni + n2 + = according to ()2ip . Propagating along eo gives 

= do{ni + 112 + 773) = -3(772 -I- 773)0- 

and thus 771 = and 772 = —773 7^ 0. Applying [Sq, ^3] to a, gives us ^n^q^, = and hence (73 = 0, 
the di derivative of which results in r2 = 0. So and qa all become zero and we have 

77o/3 = diag(0, -772,772) (n2 7^ 0) 

Qa = 0. 

This implies the vanishing of the spatial derivatives of all rotation coefficients (3q, = 0) and we 
obtain a spatially homogeneous universe of Bianchi class A, type VIq. 



^.2. n2 = na / 

First notice that a rotation about an angle (/3 exists such that n\ = — n2, ^3 = and 
"92^2 = r2 + K = 0, where 

.2 1^2 , ^2 (22) 



K = p + 3a 

This requires the existence of a solution if of the system 



doif = 

= ni + n2 

929? = sin((/?)A/ ^^ - qs 

^ dsip = cos{Lp)yJ-K - r2 



(23) 



the integrability conditions of which are identically satisfied (the first equation of (j23p guarantees 
that Oi remains zero). 

Note that the system ()23p can only give us a solution when K <Q. In that case we obtain 





2n2 





s/-K/2 


















, da — 


^/-K/2 
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and we have a LRS spacetime with 

(i) K = ^ r2 = 0: Bianchi class A, type //, 

(ii) X < =^ r2 7^ 0: when we look at the Ellis-MacCallum quantity [S] T = ^{0^ + a|), we 
find T = r| and hence this class is equivalent with Bianchi class A, type VIII. 

When K > Q (j23p has no solution, but now we can look for another rotation ip to make 
Cq, = 0. This requires the existence of a solution 99 of the system 



f 5o(/9 = 

dif = -^ + ni+n2 

d2(p = -qs 

^ dsip = -r2 



Again the integrability conditions are identically satisfied and we find 

K K , 



(24) 



Oq, = 0, Tiap = diag(2n2, 



2n2 ' 2n2 ' 



leading to the LRS Bianchi class A, type IX models. 



4.3. n2 = ns = 

These models are exceptional in the sense that they are 'purely electric' {Hat 
As in (4.2) one shows that a rotation exists under which ^3, ni and 92^2 

zero, where K is the same quantity as defined by ()22p . 
When < 0, we now obtain 











V-K/2 


















! da — 


V-K/2 




[ V-K/2 














which is a LRS spatially homogenous Bianchi class B, type VI-i model and which is 
pseudo-spherically symmetric [8]. 



0). 

r'i + K become 



When K = 0, then r2 = and we have Oq, = 0, n^^^ = 0, which is a LRS spatially homogenous 
Bianchi class A, type / model (equivalent to Bianchi class A, type VIIq [8j). 

A problem arises when K > 0: q^, ni and n2 can still be made zero by a rotation, but it is 
no longer possible to choose a frame in which all spatial derivatives identically vanish. However 
the Cartan equations can now be integrated immediately and one finds the metric 

ds^ = -dt^ + Qitfdr"^ + P{tf[de^ + sm\e)d(t>^]. 

The resulting models are the Kantowski-Sachs universes, which, as is well known, do not admit 
a 3-dimensional isometry group acting simply transitively on the t = constant hypersurfaces. 

5. Conclusion 

We have shown that the Bianchi class A dust models can be uniquely characterised as irrotational 
dust spacetimes which are purely radiative in the sense that the gravitational field satisfies 
div{H) = div{E) = under the assumption that also the magnetic part of the Weyl tensor H^/s 
is diagonal in the shear-electric eigenframe (and hence [H, E] = 0) . The only possible exception 
arises for degenerate shear and Tiaa 

= 0: when K = p + 3c7^ — > or < 0, in which case 
solutions are given respectively by the Kantowski-Sachs universes, or by the pseudo-spherically 
symmetric solutions of Bianchi class VI-i. These models however are 'purely electric' and 
cannot be called 'purely radiative' in a strict sense. 

The whole analysis presented here can be easily generalized to give a characterisation of 
spatially homogenous non-rotating perfect fluids of Bianchi class A (to be submitted). The 
next logical step is to investigate the case where [H, E] ^ Q: it is hoped that new classes of 
cosmologically interesting solutions will turn up. Alternatively one would obtain the remarkable 
result that the Bianchi class A spacetimes are the unique purely radiative ones (work still in 
progress). 
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